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Abstract 

The static and dynamic properties of the alkahne-earth atoms in their metastable state are 
computed in a configuration interaction approach with a semi-empirical model potential for the 
core. Among the properties determined are the scalar and tensor polarizabilities, the quadrupole 
moment, some of the oscillator strengths and the dispersion coefficients of the van der Waals 
interaction. A simple method for including the effect of the core on the dispersion parameters is 
described. 
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I. INTRODUCTION 



The low-lying triplet metastable states of alkaline-earth atoms have been generating in- 
creasing interest in the area of cold-atom physics for a number of reasons. One application 
is to use the ^Sq transition in calcium as a new optical frequency standard The 

use of the -^Sq ^ ^Pq transitions for fermionic ^^Sr stored in an optical lattice is expected 
to further result in an improved standard 2]. Another possible application is in the forma- 
tion of Bose-Einstein condensates (BEG) consisting of alkaline-earth atoms P, |4| in their 
metastable triplet states. The stability, size and excitation modes of BECs depends on the 
sign (and magnitude) of the scattering length, and the scattering length depends sensitively 
on the precise values of the dispersion constants 0, ^. 

Taken in conjunction, the desirability of obtaining precise values of the static and dynamic 
properties of the low-lying nsnp metastable state has greatly increased in importance. 
In this article, properties of these states are computed from valence electron configuration in- 



teraction ca, 
interaction 



culations t 
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lat use a semi-empirical model potential to describe the core-valence 
ll| . Among the data computed are the oscillator strengths for some 
of the low-lying transitions, the scalar and tensor polarizabilities, the quadrupole moments, 
and the dispersion coefficients for the van der Waals interaction between two atoms. 

II. RESULTS OF THE CALCULATIONS 
A. Methodology 

The properties of these states are computed using configuration interaction (CI) calcula- 
tions that treat the correlations between the valence particles in an ab-initio manner while us- 
ing a semi-empirical model potential to describe the core-valence interaction 
The details of this calculation are very similar to those reported in 
some minor changes in the cutoff parameters and the use of an orbital basis of larger dimen- 
sion. The polarization potentials were initially defined by tuning the potential to reproduce 
the ns, np, nd and nf binding energies of the respective singly ionized atom. The Hamil- 
tonian was then diagonalized in a basis consisting of all the two-electron basis states that 
could be formed from a set of about 140-150 single particle orbitals. The basis set contained 
orbitals up to and including ^ < 8 and the two-electron basis dimensions ranged from 1000 



apart from 



to 4000. For all practical purposes the basis for the two- valence electrons can be regarded as 
saturated. The initial binding energies obtained by this procedure were not in perfect agree- 
ment with experiment with discrepancies for the ground and excited state energies of the 
order of 0.1-2.0% (refer to ^, ^ to get an indication of the accuracy). Some further tuning 
of the cutoff parameters was done to improve the accuracy of the energy differences which 
directly impact on the accuracy of expectation values. Expectation values for multipole 
operators were computed with a modified operator that allowed for polarization corrections 



operators w 

iilQ,Q. 



The model potential is quite realistic since the direct and exchange interactions with 
the core were computed without approximation from a Hartree-Fock (HF) wave function, 
only the core polarization potential was described with a model potential. The resulting 
polarizabilities, and dispersion parameters for homo-nuclear pairs of atoms were generally 
within 0.1% of the best variational calculations for Li or Be, and for heavier atoms they were 
generally within 1-2% of results coming from large fully relativistic calculations combining 
configuration interaction and many-body perturbation theory techniques [l]| . 

The most likely source of error in the present calculations for the heavier species, Ca and 
Sr, is the neglect of relativistic effects. However, the use of a polarization potential tuned to 
the experimental binding energy will implicitly take into account the influence of relativistic 
effects upon the core electron distribution. Further, Greene and Aymar have shown that 
the spin-orbit interaction does not have maior effect on the structure of the alkaline-earth 
wave functions Li 



B. Energy levels 

The energy levels of the present calculations are given in Table H] and compared with ex- 
periment. The polarization cut-off parameters were fine-tuned to reproduce the experimental 
binding energy of the lowest states of each symmetry. In the case for states with L > 
the parameters were tuned to reproduce the center-of-gravity of the spin-orbit triplets. The 
spin-orbit splitting of the triplet states is largest for strontium and its magnitude is about 
0.001 Hartree. 

The agreement between the theoretical and experimental energy levels is sufficiently close 
to discount the possibility that energy level considerations might make a significant contri- 
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bution to the uncertainty in the oscillator strengths and polarizabilities. 



TABLE I: Theoretical and experimental energy levels (in Hartree) of some of the low-lying 
metastable states of the alkaline-earth atoms. The energies are given relative to the energy of 
the doubly ionized core. The experimental energies for the triplet states are averages with the 
usual (2 J + 1) weighting factors. The experimental data are taken from 
the 3d level for Be, Mg and Ca while for Sr it is the 4:d level. 



16 



. The md level is 



Systems 


Be 




Mg 




Ca 




Sr 




Level 


Theory 


Exp. 


Theory 


Exp. 


Theory 


Exp. 


Theory 


Exp. 




-1.011842 


-1.011850 


-0.833533 


-0.833530 


-0.660944 


-0.660932 


-0.614598 


-0.614602 


nsnp 


-0.911710 


-0.911701 


-0.733378 


-0.733788 


-0.591387 


-0.591388 


-0.547611 


-0.547612 


ns{n + l)s ^S" 


-0.774561 


-0.774552 


-0.645827 


-0.645821 


-0.517230 


-0.517228 


-0.482289 


-0.482292 




-0.739862 


-0.739855 


-0.569906 


-0.569929 


-0.485477 


-0.485478 


-0.452720 


-0.452717 


nsmd ^D'^ 


-0.729118 


-0.729113 


-0.615041 


-0.615022 


-0.568193 


-0.568180 


-0.531359 


-0.531367 



C. Quadrupole moments 

The quadrupole moment of the ^Pj" state is a static property of the state. An exact 
knowledge of its value is important since the quadrupole-quadrupole interaction has a big 
impact of the cold collision physics in metastable alkaline-earth metal atoms . Defining 
the LS coupled reduced matrix element as 

Q{L) = II J2^^C\r,) II n'n) , (1) 

i 

the quadrupole moment for a triplet state is usually defined as the moment of ^P] state 
with Mj = J; In the expression above C^(f) is the spherical tensor of rank 2. 

g = (vl/(3p;);M, = J|^r,2Co^(r,)|vl/(3p;);M, = J) , (2) 

i 

This can be written 

where the Wigner-Eckart theorem has been used twice to collapse the angular factors. The 
quadrupole moment for a level is often given for the J = 2 state. The quadrupole moment 
Q{LJ) for the J = 2 state is equal to 2{Qzz)- 
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Some older calculations of {Qzz) exist j20|, |21|. The finite element Multi Configuration 



Hartree-Fock (MCHF) calculation of Sundholm and Olsen for Be gave 4.53 au which is in 
excellent agreement with the present value of 4.54 au. The CI calculations of Ceraulo and 
Berry [2^ consistently underestimated the present quadrupole moments (e.g. 7.944 au for 
Mg) and are not listed in Table ITTl 

The quadrupole moments are compared with the recent calculations by other groups in 
Table |n] The CI+MBPT calculation 0| is a fully relativistic calculation with the post-HF 
interactions between the valence electrons and the core treated with perturbation theory 



while the interaction between the two va! 
The calculation of Santra and Greene 



ence electrons are teated with the CI ansatz. 
2^ (SG-CI) treated the two active electrons within 
a CI framework while using a model potential to represent the core-valence interaction. 
The model potential did include a spin-orbit interaction. One limitation with the SG-CI 
calculation is that it does not include the di-electronic part of the polarization potential. 

The noticeable feature of Table E] is that all three calculations agree with other with a 
total variation of less than 2%. The present results generally lie closer to the CI+MBPT 
calculation than the SG-CI calculation. The high level of agreement between three com- 
pletely independent calculations suggests that the uncertainty ascribed by Derevianko et al 
to their quadrupole moment was too big by a factor of 2. 

D. Oscillator strengths of low- lying transitions 

The oscillator strengths for the transitions to the lowest lying ^5*^, ^P^ and ^D^ states are 
given in Table |Tll The absorption oscillator strength from state ipQ is calculated according 
to the identity 

r _ 2|(V;o;Lo5' II Y.i ^»C^(fi) || jjn] LnS)\^eQn / .X 

3(2Lo + l) ■ 
The oscillator strengths for the Be triplet transitions are probably as accurate as any that 
have previously been published. The basis for the valence electrons is effectively saturated 



and the semi-empirical approach to core polarization is capable of high accuracy IjJ . For 



example, the present methodology reproduces the dipole and quadrupole polarizability of 



Be given by a close to exact calculation 23| to an accuracy of 0.2%. The present oscillator 
strengths agree very well with the experimental values given in Table |n] Not shown in the 
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Table are the ^P^ oscillator strengths of CI calculation of Weiss and the MCHF 

calculation of Jonsson et al [25] . Both of these calculations were very large and incorporated 
both core and valence excitations. The Weiss /-value was 0.447, while the Jonsson et al 
result was 0.4463. These could hardly be any closer to the present value of 0.4467. 

The present oscillator strength for the transition to the ^S*^ state in Mg, namely 0.138 is 
in excellent agreement with that obtained from the low uncertainty experiment of Andra et 
al js^, 0.139±0.003. Agreement with the large basis CI calculation of Moccia and Spizzo 
(MS-CI) is also good j27|. The MS-CI calculation is similar to the present calculation in that 
excitations are only permitted for the valence electrons. It does not allow for core-valence 
correlations so the present approach, which does, should be regarded as being more reliable. 

In the case of Ca, good agreement is achieved with the model potential calculations of 
Hansen et al 2^ for the transitions to the ^S^ and '^P^ states. A 7% discrepancy occurs for 
the transition of the 4s3(i ^D^ state. The larger difference here is expected since the 3(i orbital 
does have a tendency to penetrate into the core and therefore degrade the accuracy associated 
with model potential methods. The best ab-initio calculation is the MCHF calculation by 
Froese-Fischer and Tachiev The MCHF calculation allows for core-valence correlations 
and also includes relativistic effects using the Briet-Pauli Hamiltonian. The MCHF oscillator 
strengths listed in Table ITTl are a weighted average of the individual lines in the multiplet. 
The largest difference between the present and MCHF oscillator strengths is less than 4%. 

The multi-channel quantum defect theory (MQDT) calculations of Werji et al which 
use an i?-matrix calculation to determine the short-range parameters. Their transition rate 
data was converted to oscillator strengths using experimental energy differences and lie 
within 2-3 % of the present oscillator strengths. 



The most precise experiment for Sr is that of Andra et al j3l| which gave a lifetime of 
7.89 ± 0.05 ns for the 5p^ ^Pj^ state. This state can decay to the both the 5s5p and 5s6p 
levels and the lifetime was converted to an oscillator strength by neglecting the transition 
to the 5s6p state. This assumption is justified since the dipole matrix element will be small 
due to the (5p|6p) overlap, and the 5s6p^ ^P° — > 5p^ ^P^ energy difference of 0.0073 Hartree 
is also small. 

The comparison with the time- dependent gauge independent (TDGI) calculations of 
Merewa et al is mainly of interest because these authors also give estimates of the 
scalar and tensor polarizabilities. A quick comparison of TDGI /-values with other results 
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in Table |n] reveals that their oscillator strengths do not have the same level of accuracy 
as the other calculations. The underlying atomic structure information entering the TDGI 
formalism comes from CI calculations. 



E. The polarizabilities 

1. Theoretical treatment of polarizabilities 

This analysis is done under the premise that spin-orbit effects are small and the radial 
parts of the wave functions are the same for the states with different J. 

The Stark energy shifts for the different Lq levels in an electric field F are written as 

AE = ~aL,MoF' . (5) 

The Stark shifts for the different Mq states of the level are different and the polarizability 
is written as 

, 3Mo^ - Lo(Lo + 1) 
aioMo = ao + — Lo{2Lo - 1) ' 

where a2 is taken from the state with Mq = Lq. The total polarizability is written in terms 

of both a scalar and tensor polarizability. The scalar polarizability represents the average 

shift of the different M levels while the tensor polarizability gives the differential shift. 

In terms of second order perturbation theory, the energy shift from an electric field, F 

pointing in the z-direction is 

2 4" (^0 - i^n) 

(7) 

The polarizability can therefore be written 

V \-Mo M„y (^0 - En) 

where the Wigner-Eckart theorem has been used to isolate the M-dependent terms. Using 
the definition of the oscillator strength, eq. Q and taking the average of the energy shifts 
leads to the usual definition as a sum rule over the oscillator strengths. It is 

ao= ftLoA/o/ (2/^0 + 1) = ^ -f^ , (9) 

Mo=-Lo n On 
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where the sum includes both valence and core excitations and eon = (-E'n — En). The /-value 
distribution for the core was estimated using a semi-empirical method In this approach 
one writes 



a. 



E 



(10) 



where iVj is the number of electrons in a core orbital, is the Koopman energy, and A is an 
energy shift parameter chosen so that eq. ()10|) reproduces an accurate estimate of the core 
polarizability determined my other, independent means. 

Since the M-dependent part of the polarizability is a tensor of rank 2 and it is easiest to 
define it in terms of a2,LoLo- 



-Mo Mo 

3Mo^ - Lo(Lo + 1) 
Lo{2Lo - 1) 




(11) 
(12) 



where a2,LoLo is 

E 



1 


f Lo 


1 LA 






Lo) 



1 



3(2Lo + 1 

In terms of an /-value sum, this reduces to 

«2,Loio = " 



2|(V^0;i^0 II E^ r^C\h) II ^Pn,Ln)? 

{Eo - E^) 



(13) 



^2 9 



^n,L„=0 n,L„=l " n,L„=2 

The core does not make a contribution to the tensor polarizability since it has an equal 
impact on all the different M-levels. 

The development above is for LS coupled states, but it is common to give the tensor 
polarizability for LSJ states. These can be related to the LS states by geometric factors 
arising from the application of Racah algebra. The polarizability can be expanded 



^2 in 



]_ \ ^ fon 



(14) 



, 3M2-Jo(Jo + l) 
«JoMo — «0 H TT «2,JoJo 



(15) 



^o(^o - 1) 

where a2,JoJo is the tensor polarizability of the state with Mo = Jq. The scalar polarizability 
for the different J levels are the same and equal to the scalar polarizability in the L repre- 
sentation. The tensor polarizability between the L and J representations can be related by 
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«2,JoJo = «2,LoLo(2Jo + { ; ; ^ 4 ^ . (16) 

Lq 2 Lq 
-Lo Lo^ 

When J = 1 this reduces to a2,JoJo = — «2,LoLo/2- (This result has been checked by convert- 
ing our LS coupled /-values into LSJ coupled values and then using the standard expression 
in terms of the |(Jo || r || Jn)P matrix elements 



2. Results of calculations 

The program logic and associated numerics were initially tested by estimating the polar- 
izabilities of the ls2p level of He. The present ao of 46.6 Oq is within 0.2% of the close to 
;t calculation of Yan et al 3^. Agreement with the TDGI q;2,LoLo Rerat and Pouchan 
is not as good, but it should be noted that the TDGI calculation calculation obtains an 
cto of 49.5 Oq, indicating that the Rerat-Pouchan calculation is not quite converged. 



The present estimates of the Be polarizabilities are the most accurate that 



lave been 



3 



for ao 



published. The agreement with the Themelis and Nicolaides MCHF calculation 
is reasonable, but they give an a2,Lo-Lo that is about 20% smaller. This level of agreement 
is acceptable given that the MCHF calculation was much smaller, the 2s2p state was 
represented by a 3 configuration MCHF wave function while 14 configurations were used to 
represent the excited states. 

Only a moderate level of agreement is achieved with the TDGI polarizabilities for Be, 
Mg and Ca HHQ- 

The static polarizabilities agree at the 10% level while the TDGI 
estimates of the tensor polarizability are up to 50% different. The lower level of accuracy 
achieved by the TDGI calculations is consistent with the earlier discussion concerning the 
accuracy of the oscillator strengths. 

A recent measurement of the tensor polarizability for the state of Ca using an atomic 
polarization interferometer gave 2.623 ± 0.015 kHz/(kV/cm)^ or 10.54 ± 0.06 Oq 0. The 
tensor polarizability of the J = 1 state is determined from the Ca entry in Table |H] by 
multiplying by -i according to eq. (fTBj) . The present calculation gives a2,JoJo = 14.2 aj] for 
the 3Pi° state. A very early estimate of the tensor polarizability for this state was 12.9 ± 3.2 
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Uq j5l| and another independent experiment gave a2,JoJo = 12.1 ± 0.8 Og [5^ 

The scaler polarizabihty of the state has not been measured directly, but there have 
been measurements of the difference between the polarizabilities of the 4s^ ^S^ ground state 
and the {m = 0) ground state. Morinaga et al ^| obtained 90.4 ± 13.5 ajj for the 
difference in the polarizabilities. Using the polarizabihty of 159.4 for the Ca ground state 



Ul, and the present '^P^ (m = 0) polarizabihty of 295.3 - 2 x 14.2 = 266.9 al gives 107.5 
for the difference in the polarizabihty. 

The Stark frequency shift of Li and van Wijngaarden of 12.314 ± 0.041 kHz/(kV/cm)^. 
for the 4s^ ^Sq ^P" (m = 0) transition |5J| converts to a polarizabihty difference of 
98.98 ± 0.33 al 

Taken together, present estimates of a2,JoJo larger than experiment by about 20% 
while estimates of the ao(4s^ ^S*^) — ao(4s4p ^P") polarizabihty difference are about 10% 
too large. Rectifying the situation in a non-relativistic calculation could be problematic 
since an improvement in a2,jo,Jo will result in the theoretical polarizabihty difference drifting 
further away from the experimental polarizabihty difference. 

The obvious improvement that could eliminate this problem would be the inclusion of 
the spin-orbit interaction. The largest contribution to the polarizabihty comes from the 
transitions to the levels. The spin-orbit splitting leads to the excitation energies for 
4s3(i states with differing J fluctuating by about ±2%. Given the cancellations that occur 
in the evaluation of eq. (|T3|l it is possible that introduction of spin-orbit splitting could lead 
to a Ca tensor polarizabihty in better agreement with experiment. 

There have been no measurements of the tensor polarizabihty for the other alkaline-earth 
atoms. This should be rectified since it would be a very useful diagnostic with which to 
assess the accuracy of the structure models of the metastable states. 



3. Alternate treatment of core 

It is desirable to partition the core /-value some into contributions that arise from exci- 
tations to final states with different core-|-valence angular momentum, L^- Therefore, it is 
possible to write symbolically 

encore ^ ^ Oicore,Lx ■> 



L 



T 
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where acore,LT wiU include all the contributions from the different magnetic sub-levels, i.e. 

(18) 

Mt 

For any of the core dipole excited magnetic sub-levels one can write 

_^ f{00:LM^n;LTMT) 



'On 



The final states, \LtMt) can be expanded in terms of uncoupled states, e.g. 

\LtMt) = ^{lmLM\LTMT)\lmLM) , (20) 

mM 

where Im refers to the angular momentum of the excited core and LM refers to the angular 
momentum of the metastable state which is acting as a spectator. Therefore, it is 
possible to decompose the oscillator strength as 

/(0;OOLM^n;LTMT) = ^|(lmLM|LTMT)|LTMT)|V(0;OOLM^n;lmLM) . (21) 

mM 

The polarizability can also be expanded in terms of uncoupled states 

J2 I {lmLM\LTMT)\LTMT)\^acore,mM/{2L + 1) . (22) 
mM 

The factor of (2L+1) in the denominator arises due to the sum over spectator states. We now 
assume that the excitations for the core occur independently of the state of valence electrons 
which act as spectators. Therefore, the contribution to the polarizability is independent of 
M. Further, the core initially has a net angular momentum of zero and therefore the 
different magnetic sub-levels of the core excitations should give equal contributions to the 
polarizability, hence 

_ Oicore /rjq\ 
(^core,mM — ^ ■ {^"^J 



The final result is 

|(lmLM |LrMr) I-LtMt) r — 

mM 

which can be simplified by summing the Clebsch-Gordan coefficients to give 



aLrMr = J] | (lmLM|L^M^) |LtMt) T-^^ , (24) 



_ (2Lr + l)acore /r,r\ 
acore,Lr- 3(3^+1) " ^^^^ 
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When particular values are substituted into eq. the distribution of the core /-value 
sum into the ^S°, and manifolds is given in the proportion | : | : |. This is of 
course just the statistical weighting associated with the {2Lt + 1) degeneracy factor. It is 
simple to verify that such a proportion means the net contribution of the core to the tensor 
polarizability as defined by eq. p4|l is zero. 



F. The van der Waals coefficients 

The van der Waals coefficients given in Table |n] are those for a pair of ^Pq states. The 
dispersion parameter, Cq is simple to compute since both of the atoms have a net angular 
momentum of zero. The expression is 

Ce = - V fo,nJo,n, ^ _ ^26) 

2 eomeoriaieoni + eona) 

The present dispersion parameters are slightly larger than those of the SG-CI calculation of 
Santra and Greene [22]. Taking the case of Sr, the difference here is about 15%. About half 
of this difference can be attributed to the core since Cq = 5668 au when core excitations are 
omitted from eq. ()26|) . So part of the discrepancy arises from the neglect of the core in the 
SG-CI calculation. One cautionary note should be made. Santra and Greene reported Cq for 
the 3Po° state. Since the ^Pq state is the most tightly bound state of the 5s5p multiplet one 
expects the present LS coupled calculation to have a slightly larger Cq. The quantitative 
impact of spin-orbit splitting can best be determined by separate evaluations of Cq for the 
J = 0, 1 and 2 states. 

The van der Waals coefficients that are relevant to BEG studies are those between two 
states. The algebra related to this is somewhat messy and the coefficients are presented 



in the formalism of Santra and Greene 
between two ^P] states is defined as 







22]. The intermediate dispersion coefficient 



R _ (_'\\Jl-J2+l /o,ni/o,n2 fr,^\ 



ni,n2 



where rii has angular momentum Ji and n2 has angular momentum J2. 

This LS coupled oscillator strengths were converted into the LSJ coupling scheme using 
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the identity 

fiJo^Jn) = fiLo^Ln)i2Lo + l){2Jn + l){ "} • (28) 




When the sum, eq. (j2ZI), was evaluated, the core /-value distribution was included using 
eq. (j25j) to partition it into ^5**^, ^P'^ and ^D'^ excitations. 

The results of our calculations are presented in Table IIIII and compared with earlier 
CI+MBPT calculations of Derevianko et al and the SG-CI calculations j^]. There is no 
apparent experimental activity on the metastable states of Be and the present data in the 
Table were only included for reasons of completeness. 

The present calculation and the CI+MBPT calculation could hardly be in any better 
agreement for magnesium. The largest disagreement for any of the Bj^ j^ coefficients was 
1.2% for the i?2,2 coefficient. Agreement with the SG-CI calculation is not as good with the 
occasional discrepancy of 5% and it is noticeable that the present and CI+MBPT results do 
tend to be larger in magnitude. 

For calcium there is a tendency for the present results to be from 1% to 5% larger in 
magnitude with the differences being smaller for the larger values of Ji and J2. The present 
dispersion coefficients all lie within the 10% uncertainty that Derevianko et al associate 
with their results. Agreement with the SG-CI calculations is not so good with discrepancies 
exceeding 10% being common. 

The pattern for strontium is similar to that seen for calcium. The present Bj^^j^ coef- 
ficients are larger than the CI+MBPT data for Bi i and smaller for -63,3. The differences 
with the SG-CI calculation are generally larger than those with the CI+MBPT calculation. 

Some general trends are noticeable. The SG-CI calculation always gave the smallest 
result for Bi i, -82,17 B31, -83 2 and -83,3. Furthermore, the sum Bj- j. for the present 

calculations and CI+MBPT calculations are consistently bigger than the SG-CI calculations, 
with the difference becoming larger as the atom gets heavier. This could be a manifestation 
of the increasing importance of the core contribution to the B j.^j. coefficients as the atom 
gets heavier. 

It is also evident that some of the uncertainty estimates of the SG-CI calculation were 
somewhat optimistic. For example, they give Bi i = 139 + 7 au for strontium. The contribu- 
tion of the core /-value sum to this dispersion parameter is 12.0 au. So the core contribution. 
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which is not incorporated in the SG-CI calculation, is larger than their estimated uncertainty. 



III. CONCLUSIONS 

A systematic study of the properties of the alkaline-earth atoms reveals that the present 
non-relativistic approach reproduces the results of the CI+MBPT ansatz of Derevianko et 
al yi to better than 5%. Agreement with the model potential CI calculation of Santra 

ri 

and Greene j22| is not so good with discrepancies of 10-15% occuring for the spherical 
part of the Cg dispersion coefficient. Due to the unknown impact of the spin-orbit energy, 
splitting upon the polarizabilities and dispersion coefficients, it is not possible to make a 
definitive statement about any reasons for the differing levels of agreement with these two 
other calculations. However, we do suspect that the omissions of the SG-CI model, i.e. the 
di-electronic two body polarization potential, the non-usage of a dressed dipole transition 
operator, and the lack of core excitation terms in the dispersion sum rules all contribute in 
part to the differences with the SG-CI model. 

It should be noted that previous studies with the present model for the alkali atoms and 
singlet states of the alkaline atoms demonstrated that the method could predict a number of 



expectation values with an overall accuracy of 1-2% or better 1^. The presence of spin-orbit 
energy splitting, and the existence of a ^D'^ state very close in energy to the metastable 
level leads to a decrease in accuracy for atomic properties such as the tensor polarizability 
that are sensitive to these energy differences. Additional high precision measurements of 
the tensor polarizabilities for the Mg, Ca and Sr would be certainly be worthwhile since 
the sensitivity of this parameter to the fine details of the wave function should help in the 
refinement of the two-body potentials used to characterize ultra-cold collisions. 

It is interesting to speculate whether the better agreement with the CI-I-MBPT cal- 
culations could be achieved by incorporating a spin-orbit potential into the Hamiltonian 
and using jj coupling. Alternatively, a fully relativistic treatment, using a relativistic HF 
wave function might be necessary. Resolution of these questions requires that explicit cal- 
culations be made to determine the additional physics needed to eliminate the anomalies 
between the present calculations and experiment and between the present calculations and 
the CI-FMBPT calculations. 
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TABLE II: Properties of the metastable levels of the alkaline-earth atoms and He (note, the 
lowest "^P" level is not metastable in He). The oscillator strengths to the lowest ^S^, ^P^ and 
^D^ states are given as fi^L^). The scalar and tensor dipole polarizabilities are ckq s-nd a2,Loio 
respectively. The quadrupole moment Q is given for the ^P2 state while the dispersion parameter 
Cg is that for two ^Pq states. The He "Other Theory" row reports the results of close to exact 
calculations with the exception of a2,Loio- "^^^ present oscillator strength to the He ^D^ state is 
not to a physical state, rather it is to the lowest energy pseudo-state. All quantities are in atomic 



units and the numbers in brackets are the uncertainties in the last digits. 



Met no a 


J{ ^ ) 




/(, ^ ) 


Q 


ao 












He 










Present 


-0.1797 




0.6251 


10.264 


46.66 


69.62 


5102 


Other Theory 


-U. i ( y ( \oo\ 




n ^^1 no FqqI 
O.DlUz ^ol 


1 n ofip; Fq/iI 
lU.ZDO 


40. r 1 


o/.uy ^pj 










Be 










Present 


0.08187 


0.4467 


0.2948 


4.54 


39.02 


0.558 


220.3 


MCHF [37] 










39.33 


0.47 




TDGI 


0.026 




0.154 




36.08 


1.04 




B-spline CI [391 


0.0823 


0.453 


0.295 










Experiment 


0.089(3) [4^1 


0.44(2) [^ 


0.29(1)[4J^ 
















Mg 
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0.1383 


0.6167 


0.6287 
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101.9 


-14.24 


1004 










8.46(8) 








SG-CT [221 








8.38 






9S0f30~) 


TDGI [42] 
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0.625 




90.7 


-19.64 




IVIS-CI [27] 


0.1354 


0.6383 


0.6336 










Experimental 


0.139(3) [26] 


0.55(4) \M 


0.62(4) [44] 
















Ca 










Present 


0.1582 


0.5071 


0.08136 


12.96 


295.3 


-28.36 


3363 


CI+MBPT [4] 








12.9(4) 








SG-CI [22] 








12.7 
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TDGI [4^ 


0.163 
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-50.0 




Cl+model [2^] 
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0.0873 










IVICHF [29] 
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0.0806 










Experimental 


0.12(2) [46] 
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Sr 
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IvIQDT [22] 
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Experimental 
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TABLE III: The intermediate dispersion coefficients, Bj-^ j^ for two alkaline-earth- metal atoms 
in the metastable state. The ^ \Bj.^j.\ column sums the absolute value of all the entries in 
each row (with off-diagonal elements added twice) . The numbers in brackets after the data are the 



uncertainties ascribed to the CH-MBPT and SG-CI calculations. 



Method 


Bi.i 


B2,l 


B2,2 






Bz,2 


^3,3 


Eis./.,./, 


Present 


-6.901 


11.29 


-27.77 


Be 


-13.52 


22.31 


-27.21 


220.3 


Present 
CI+MBPT [4] 
SG-CI |22J 


-37.46 
-37.8(38) 
-35.6(2) 


43.90 
43.9(44) 
42.5(2) 


-52.73 
-52.1(52) 
-51.9(2) 


Mg 


-76.49 
-76.7(77) 
-73.4(7) 


90.76 
90.1(90) 
88.6(7) 


-157.2 
-156.4(156) 
-152(2) 
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1002 
976 
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CI+MBPT [4] 
SG-CI [22] 


-96.95 
-91.7(92) 
-81(3) 


130.7 
123(12) 
119(5) 


-176.8 
-167(17) 
-176(8) 


Ca 


-233.5 
-225(23) 
-203(10) 


317.5 

306(31) 
302(20) 


-604.8 
-600(60) 
-553(70) 


3363 
3250 
3087 


Present 
CI-I-MBPT [4J 
SG-CI [22j 


-165.9 

-158(16) 
-139(7) 


213.9 

203(20) 
196(9) 


-278.9 
-264(26) 
-280(10) 


Sr 


-416.9 
-415(42) 
-370(30) 


556.0 

555(56) 
546(50) 


-1231 
-1290(130) 
-1210(200) 


6074 
6090 
5780 
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